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Abstract 

The foUowing aspects of hadronic final states in deep inelastic lepton scattering are reviewed: 
measuring as from multi-jet production rates and event shapes; alternative jet algorithms for DIS; 
power-suppressed corrections to event shapes; comparing jet fragmentation in e'^e~ annihilation 
and DIS; final states in the BFKL and CCFM formulations of small- 2; dynamics; exotic (instanton- 
induced) final states. 



Introduction 



This talk will review a selection of topics concerning 
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hadronic final states in deep inelastic lepton scattering 
(DIS). One of the primary aims of studies of DIS 
final states is to test the predictions of QCD in 
I more detail than is provided by measurements of the 
^otally inclusive structure functions. By analogy with 
(-H 'e~^e~ annihilation, it is natural to begin by studying 
• • various global measures of the jet-like properties of 
.5_J the final state. This is the topic of the following 
section. Amongst the quantities considered are the 
5-H multi-jet production rates, defined according to various 
. 5^ jet algorithms, and event shape parameter distributions. 
Both of these have to be defined in an infrared- 
safe way, so as to be insensitive to the long-distance 
behaviour of QCD and hence calculable in perturbation 
theory. Recently, it has been found that observables 
satisfying these conditions can still diiTer substantially 
in their sensitivity to non-perturbative physics, through 
the magnitude and energy dependence of corrections 
that are suppressed by inverse powers of a large 
momentum scale. Some recent theoretical ideas on 
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this question and their predictions for DIS will be 
discussed. Experimental and theoretical results on the 
fragmentation spectra of DIS current jets will also be 
reviewed. 

Section 3 deals with final state features of DIS in the 
region of very small Bjorken x values currently being 
investigated at HERA. Here the theoretical analysis is 
more difficult and a number of basic points remain to 
be clarified. Some of the most promising indicators of 
possible new dynamics at small x concern final state 
properties, such as the transverse energy flow and the 
production of forward jets with transverse momenta 
comparable to the DIS momentum transfer Q. 

In section 4, a more speculative but exciting 
possibility for DIS final states is discussed, namely non- 
pcrturbative processes that might be induced by QCD 
instantons. The associated final states would have high 
hadron multiplicities without any prominent multi-jet 
production. 

Section 5 contains some concluding remarks. 

2. Jet Physics 

2.1. Jet rates 

Multi-jet production rates are at present the only 
features of DIS final states that have been fully predicted 
to next-to-leading order in QCD perturbation theory . 
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Figure 1. Typical diagrams for neutral current deep-inelastic 
scattering: (a) Born term, (b) QCD Compton scattering, (c) 
boson-gluon fusion. 



Figure 2. Results from the ZEUS collaboration on as for 
different regions. 



Denoting by cr„+i(x, Q^; j/cut) the cross section for the 
production of n jets plus the remnant jet, at given values 
of Bjorken x, momentum transfer-squared and jet 
resolution j/cut, we have schematically 

+ <(g)B„(x,g2;y^^^) + ...(i) 

where An and _B„ are calculable in terms of the parton 
distribution functions of the proton. The leading terms 
An are given by the corresponding tree graphs as 
illustrated in figure^ Once the next-to-leading term i?„ 
is known for n > 1 , this provides a method for measuring 
the strong coupling ag. 

There are several different ways of defining jets, with 
correspondingly different definitions of the jet resolution 
2/cut- So far, the ncxt-to- leading function B2 for 2+f 
jets has been computed only for the modified JADE]et 
clustering algorithm |l|, ||. Here one defines for each 
pair of particles or clusters i and j the quantity 

= 2E,E.j (1 - cos % ~ ml /W^ (2) 

where rriij is the invariant mass of the pair and W 
is the overall hadronic centre-of-mass energy {W'^ — 
(5^(1 — x)/x). The pair with the smallest value of Uij 
are clustered, the process being repeated until all j/^ 's 
are above ymt- The proton remnant is included in the 
clustering procedure as a 'pseudoparticle' carrying any 
missing longitudinal momentum that is lost down the 
beam pipe. After clustering, the clusters remaining at 



resolution ycut are defined as jets, the one containing the 
pseudoparticle being the remnant jet. 

Using the above jet definition at j/cut — 0.02, 
the HERA experiments have obtained evidence of the 
running of as , and have measured its value as ^ 

as(Mz) = 0.f23±G.0f8 (Hf), (3) 

as(Mz) = 0.ff7±0.GG5i°:°^^± 0.007 (ZEUS), (4) 

where the three components of the error in are 
statistical, experimental systematic and theoretical 
systematic, respectively. In figure ^, for example, the 
ZEUS results on as as a function of are shown, 
together with the expected behaviour for various values 
of the 5-flavour MS QCD scale A^. The inner error 
bars show the statistical errors only, while the outer 
ones represent statistical and systematic errors added 
in quadrature. 

In addition to the modified JADE algorithm used 
in the above analyses, there are other possible jet 
definitions, which have some theoretical advantages but 
are not yet predicted to 0{a^). For future reference I 
mention here the factorizable JADE algorithm with 

= 2E,Ej{l - cos e,j)/Q^ ~ ^ 

i.e. differing from (H) only in the normalization, and the 
k±-algorithm for DIS Q, in which 

y,, = 2min{i?f , E^}{1 - cos^^O/Q^ ~ kl,^/Q^ , (6) 

where k± represent the transverse momentum of i 
relative to j, or vice versa, whichever is the smaller. In 



2 



both these algorithms the definition of the remnant jet is 
also different from that in the modified JADE algorithm: 
instead of the pseudoparticlc, one introduces a fixed 
momentum vector pr — xp, where p is the incoming 
proton momentum, and computes ytr, obtained by 
clustering pi with pr, along with j/y at each stage. If 
the smallest of all these is ykr, then k is classified as 
part of the remnant jet and is not available for further 
clustering. 

A substantial part of the theoretical systematic error 
in the as determinations (||) and (|^) is associated with 
hadronization corrections, i.e. the non-perturbative 
corrections applied to the perturbative predictions 
before comparing them with hadron-level data. The 
/e_L-algorithm in particular may have some advantages 
in reducing this source of error. As we shall see in 
section 2.3, there are theoretical indications that non- 
perturbative effects should be smaller for algorithms 
based on the /c^-resolution variable rather than the 
JADE resolution (|) or (|). 

2.2. Event shapes 

Shape variables that describe the "jettiness" of the final 
state have proved useful in e"'"e~ annihilation studies; 
one of the earliest and still most commonly-used is the 
thrust, T, defined as 



T = max 



(7) 



where the sum is over all final-state particles and the 
maximum is with respect to the direction of the unit 
vector n. The thrust has value one-half for a fully 
isotropic final state, and its value approaches unity as 
the configuration in the hadronic centre-of-mass frame 
becomes more two-jet-like. 

If the definition (0) is taken over directly to describe 
DIS, the value of the thrust is strongly affected by the 
properties of the proton remnant. This is not very 
satisfactory since the remnant is not involved in the 
hard scattering and in any case is often not seen in the 
detector. An alternative definition involves looking at 
the final state in the Breit frame of reference instead of 
the hadronic cm. frame. The Breit frame is the one in 
which the momentum transfer from the lepton is purely 
spacelike and lies along the z-axis: say, = (0, 0, 0, Q). 
In this frame the current jet is usually in the same 
hemisphere as (the current hemisphere), while the 
remnant jet is in the opposite (remnant) hemisphere. 
We may therefore define the "current jet thrust" T^, as 
2P/Q where P is the total longitudinal momentum in 
the current hemisphere in the Breit frame. 

The distribution of the current jet thrust has not yet 
been calculated to next-to-leading order. A calculation 
and experimental measurements of this quantity would 



be valuable for as determination and for comparison 
with e+e^ results. We shall see in the following 
subsection that recent ideas about non-perturbative 
contributions could be tested by such a comparison. 

One shape variable whose distribution can already 
be predicted to next-to- leading order is 1/2+1, the value 
of the jet resolution ^cut at which two jets plus the 
remnant jet can just be resolved in the final state. This 
distribution is essentially just the derivative of the 1-1-1- 
jet rate with respect to j/cut, which can be deduced (for 
the modified JADE jet algorithm) from the jet rate 
calculations described above. So far, the data have 
not been presented in this form. Again, we shall see 
below that ideas about non-perturbative contributions 
could be tested using this quantity, by studying the Q^- 
dependence of the discrepancy between the perturbative 
prediction and experiment. 

2.3. l/Q corrections 

A field of recent theoretical activity that needs more 
experimental input, which could be provided by studies 
of final states in DIS, involves the investigation of l/Q 
corrections to hadronic event shapes ^, ^ |ll|, ^ 

In e"^e annihilation, it is well known that 
many event shape variables receive significant non- 
perturbative contributions of the form X/Q, where 
A is typically of the order of 1 GeV and Q is the 
hard process scale, which in e~^e~ annihilation is the 
centre-of-mass energy, Q = Ecm- This is seen for 
example in the mean value of the thrust (figure 
the discrepancy between the data and the perturbative 
prediction (dashed) shows a clear 1/Ecm dependence. 
As also illustrated in figure H, the full dependence 
on Ecm is well described by the Monte Carlo event 
generators JETSET and HERWIG The 
discrepancy has therefore customarily been described as 
a "hadronization correction" , estimated according the 
the models of the hadronization process that are built 
into those programs. Similar corrections, also with a 
l/Q dependence, are found in differential event shape 
distributions. 

Recent theoretical studies suggest that l/Q correc- 
tions are not necessarily related to hadronization, but 
may instead be a universal soft gluon phenomenon as- 
sociated with the behaviour of the running coupling at 
small momentum scales [^l[ The term 'universal' 
means that they could be expressible in terms of a few 
non-perturbative parameters that are not themselves 
calculable but have calculable process-dependent coef- 
ficients. 

Final states in deep inelastic scattering are an 
excellent place to test these ideas experimentally. If 
the conjectured universality is true, there should be l/Q 
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Figure 3. Mean value of 1 — T in e+e annihilation. 



corrections to shape variables for DIS final states, where 
Q is now the usual DIS momentum transfer variable. 
Furthermore the coefficients of 1/Q should be related to 
those measured in e+e~ annihilation p^ . 

From the experimental point of view, DIS appears 
to be a good process in which to study 1/Q effects, 
since the wide range of values available in a 
single experiment should make it straightforward to 
disentangle the relevant power-behaved terms from the 
logarithmically-varying perturbative terms. Note also 
that the terms in question are dominant over the 
more familiar higher-twist corrections, which behave like 

Before giving some illustrative predictions for DIS 
final states, let us recall the mechanism proposed 
in ref. as the source of 1/Q corrections. For 

definiteness, consider again the mean value of the thrust 
in e~^e~ annihilation. The "improved" leading-order 
perturbative prediction of this quantity is of the form 



(i-r> 



Q 



dkj_asikj_)fTikj_,Q) 



(8) 



By "improved" we mean that in the perturbative 
prediction we use the running coupling constant 
evaluated at a scale given by the transverse momentum 
fcj^ of the emitted final-state gluon [|l9|. The function 
fx turns out to have the behaviour 

fTik^,Q) ^ iCp/nQ for fc^ « Q . (9) 

Substituting this and the perturbative expansion for 



Figure 4. Solid curve: possible infrared-finite behaviour of the 
strong coupling at low scales. Dashed: one-loop running as{k±). 
Dot-dashed: expansion of one-loop as{kj_) to second order in 
as {Q) for Q/A = 100. 



ar*(fc±) = as ((?) + & In 



Q 



HQ) 



{b = [33 - 2A^/]/67r) into eq. (§, we find 
4Cf 



(l-T) 



-as(Q)^n! [basiQ)Y 



(10) 



(11) 



Because of the n! in the coefficient, the series is strongly 
divergent. Therefore by "improving" the perturbative 
prediction we have made it meaningless! 

The reason for the divergence is that the perturba- 
tive expression for as{k±) diverges at the Landau pole, 
k± = A. Therefore, if we use as™* in 6q. (^), the integral 
is not well defined. On the other hand if we truncate 



the series (10) at any finite n there is only an integrable 
divergence at k± = 0. Consequently, for consistency, 
the series must diverge. 

One can attempt to give a meaning to the integral 
involving as°'* by specifying how to deal with the 
Landau pole. The various options differ by terms 
of order 1/Q and therefore one is led to expect a 
term of this form in the true answer, whose coefficient 
cannot be determined within perturbation theory. This 
is an example of the (infrared) renormalon ambiguity 

In ref. Oil a more phenomenological approach is 
proposed. Suppose that a meaningful, universal form of 
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CKs(fc±) can be defined for all values of k±; for example, 
the solid curve in figure |[ Then the integral 



(12) 



exists for all > 0. For /z sufficiently large :s> A) we 
have asik±) — o:s^^^{kj_) for k± > /x. Therefore we may 
write 

(l-T) ~ (1-T}P"* + (1-T)P°" (13) 

where (1 — T)^'^'^^ represents the second-order perturba- 
tive prediction, while (for A ^ ^ ^ Q) 



(1-T) 



pow 



4Cf 



dk± 



as{k_i_) 



-as{Q) - bin 



Q 



HQ) 



nQ 



ao{fi) - as{Q) 



M In^ + 1 aiiQ) 



(14) 



The negative terms in the integrand express the fact 
that we have subtracted the perturbative contribution 
from k± < fj, and replaced it by the non-perturbative 
expression (p^. As shown in figure ^, the subtracted 
part is probably smaller than ag{kj_) and therefore we 
expect a positive power correction. Good agreement 
with the thrust data shown in figure ^ is found for 
ao(2 GeV) = 0.52 ± 0.04 |l|. 

Using the same approach, one finds power correc- 
tions like ( |l^ for a wide range of final-state observ- 
ables in hard processes, including deep inelastic scatter- 
ing [15l Era . In general the correction is of the formm 



4C. /mV, . 
a 7^ log 



-as{Q)-b{\n^ + -]aliQ) 
/i r 



(15) 



where a, r and s (s = or 1) are constants that depend 
on the observable F, Ci is the relevant colour factor 
{Ci = Cp = 4/3 for the quantities we are considering) 
and 



dkj_k^j^ ^as{kj_) 



(16) 



Predictions for the coefficient a and the exponents r 
and s for various e"'"e~, deep inelastic and Drell-Yan 
observables are listed in table 1 p^ . 

f In ref. there is a constant added to the final 1/r, 

corresponding to the use of a difl'erent renormalization scheme 
for the definition of a^-i. 



F 


a 


r 


s 


\ 1 CC 


-1 


1 







37r/2 


1 





(0"L/0-tot)0C 


7r/4 


1 







1 


1 





(S>oc 


1 


1 


1 




1 


1 







- 


2 


- 




-1 


1 







/ ^ 
V i-i 


1 





(S/2+l>DIS 


1 


1 





(S/2+l>Dis 


- 


2 


- 


('3i/Q>DY 


1 


1 


1 






2 





Table 1. Predictions for power corrections to e'^e , deep 
inelastic and Drell- Yan observables. 



The e^e~ quantities T, C and (Tl are the 
thrust, C-parameter and longitudinal cross section, 
discussed in ref. [|ll|. Mj represents the jet mass- 
squared. In one-loop order this may correspond to 
the heavy jet mass-squared or the heavy-light jet 



Mf , there being no difference between 



difference 

these quantities until multi-jet contributions become 
significant. Discrepancies between the coefficients of 
1/Q corrections to these quantities therefore measure 
the non-perturbative effect of multi-jet contributions not 
associated with the running of as [ |l3| . In our treatment 
the effective expansion parameter here is ao ~ 0.5, and 
so the predicted coefficients would not be expected to be 
more reliable than about ±50%. Similarly for J5, which 
represents either the total or wide jet broadening |^ at 
one-loop order. 

Quantities like {B)^^ have enhanced (s = 1) leading 
power corrections associated with the logarithmically 
growing phase space for gluon emission with limited 
fcj^ . Such quantities also have non-enhanced terms with 
the same power, but these are more difficult to predict 
because they receive contributions from the non-soft, 
wide-angle region. 

The quantities {ysYl^ and (ys)'^^ are the values of 
the jet resolution j/cut at which three jets are just 
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resolvable, using the JADE and fc^ -algorithms. The k±- 
algorithm is "better" in that its correction is suppressed 
by one extra power of Q. This could be related 
to the observation that this algorithm has smaller 
hadronization corrections. 

Note that predictions for the coefRcient a and 
the logarithmic exponent s are not shown in table 1 
for observables which, hke (ys)^^^, have no 1/Q 
correction. When the leading power correction is of 
order there are other mechanisms that can give 

contributions of the same order, making the predictions 

[1. 

observables for deep inelastic 



less straightforward 1 14 
Turning 



to the 

scattering, we consider first the current jet thrust 
as defined previously. For this quantity the power 
correction (unlike the perturbative contribution) is 
predicted to be approximately equal to that for the 
thrust in e^e^ annihilation at -Ecm = Q, independent 
of Bjorken x. As discussed above, corrections to this 
relationship are higher-order in ag, and so the equality 
should be good to about ±50%. 

The observable (7/2+1)013 mean jet resolution 

at which two jets plus the remnant jet are just resolvable 
in DIS using jet clustering algorithm A. We denote by 
A ~ mJ, fj and k± the modified and factorizable JADE 
and k± algorithms defined above. The mJ algorithm is 
the only one for which next-to-leading order predictions 
are available at present Note that in this case the 
factor of ^ xjiX — x) in the coefRcient a means that the 
power correction should be approximately the same as 

W ^Q^{\~x)lx, 



that in e annihilation at 



rather than at £^cm = Q- 

In the factorizable JADE-type algorithm (fJ), on the 
other hand, the scale of the power correction is Q instead 
of W ^ making it larger than that for the modified JADE 
algorithm. One should remember, however, that this 
is also true of the perturbative terms, so the relative 
correction is similar. 

As in e^e~ annihilation, the /cj^-type algorithm for 
DIS 1^ has a smaller (order 1/Q^) correction, suggesting 
a smaller "hadronization" correction. 

For Drell-Yan processes, qt in table 1 is the 
transverse momentum of the lepton pair and Q is its 
invariant mass. The mean value of qtlQ is expected 
to have a 1/Q correction, again log-enhanced due to 
the growing phase space for gluon emission. The 
possibility of a XjQ power correction to the Drell-Yan 
cross section itself has also been considered but seems 
less likely @ [|. 

2.^. Fragmentation studies 

In addition to jet rates and event shapes, one can 
investigate the fragmentation properties of jets in DIS 
final states p3, 031. Here again it is advantageous 
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Figure 5. Mean charged multiplicity in e+e (open symbols) 
and DIS final states. 



to study the current hemisphere in the Breit frame, 
which should be comparable to a single hemisphere 



of an 



final state at Ecm = Q- In particular 



the average multiplicity of charged hadrons should be 
asymptotically independent of Bjorken x, with a Q^- 
dependence given by the next-to-leading-logarithmic 
(NLL) prediction M 



(rich) Aas(Q) exp[c/Va^(Q)] , 



(17) 



where the normalization A is not predicted (but should 
be equal to that for a single e~^e~ hemisphere) while 
the exponents b and c are known constants. As shown 
in figure |[ the data on charged multiplicity in current 
fragmentation do follow the Q-dependence of the e+e~ 
data, which in fact agree well with the QCD prediction 

(0). 

One can go on to compare other properties of 
the charged multiplicity distribution in current jet 
fragmentation with the corresponding results for a single 
hemisphere in e+e^ annihilation. For example, the 
second binomial moment {nch{nch ~ 1)) should satisfy 
the x-independent relation p5|| 



{nchjnch - 1)) 
{rich) 



55 5 

Nf 

8 81 ^ 



A 



as(Q) 



+0[a,{Q)] 



(18) 



The data show that the shape of the distribution does 
appear to depend on and not explicitly on x, but a 



comparison with eq. (18) has not yet been presented. 
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Figure 6. Position of the peak in the ln{l/xp) distribution for 
e+e~ (open symbols) and DIS final states. 



Studies of the current jet fragmentation function, 
i.e. the inclusive single particle momentum distribution, 
are also being performed ^ . One expects that the 
distribution of the logarithm of the momentum fraction, 
^ = ln(l/a;p), where Xp = 2p/Q in the Breit frame 
current hemisphere, should be the same as that in e+e~ 
annihilation at i?cm — Q- The distribution in ^ should 
be approximately Gaussian, with a peak that moves 
linearly with InQ. More precisely, to NLL accuracy we 
expect the peak to be at ^ = ^* where 



1 



■In 



A, 



off 



6tt 



(19) 

Here A^g is an effective scale (not necessarily equal to 
A^) and the ellipsis represents C'[q;s(Q)] corrections. 
As shown in figure |^, a linear leading In Q dependence 
is seen, in good agreement with the e'^e" data and in 
disagreement with a simple phase space model. 

Again, it will be interesting to test other QCD 
predictions concerning the form of the ^-distribution. 
For example, the NLL result for the r.m.s. width a is |2q] 



36 



Q_ 

Acff 



/ 6n 33 1 
as(Q) 8 +4 ^'^ 



(20) 

As the data become more precise, it will be necessary 
to go beyond the NLL treatment, and systematic 
differences between e~^e~ and DIS jet fragmentation 
should become manifest, due to the different higher- 



order corrections, quark flavour composition and 
threshold effects in the two processes. 

3. Small X Final States 

3.1. Theoretical remarks 

One of the areas of greatest current interest in DIS 
physics is QCD dynamics in the region of small 
Bjorken x. Although HERA has opened up a much 
larger window on the small- a; region, it is still the 
case at present that different theoretical approaches 
can account for the increase in the totally inclusive 
structure function F2 at small x It is therefore 

important to consider less inclusive aspects of DIS, 
and in particular to identify features of the final state 
that might distinguish between the different mechanism 
proposed for the small- a; rise in F2. 

Outside the small- a; region, the standard approach 
embodied in the DGLAP Q^-evolution equations [p8| 
has a firm theoretical basis in the operator product 
expansion, renormalization group and factorization 
theorems. A simple extrapolation of DGLAP evolution 
to small X still appears to be consistent with the 
structure function data. The corresponding final states 
are dominated by strongly ordered configurations of 
transverse momentum kt and virtuality k^ along the 
ladder of partons connecting the soft, low- virtuality 
proton constituents to the hard subprocess (figure 
Such configurations generate the greatest number of 
logarithms of in the evolution of the structure 
functions. Sub-leading corrections are classified 
according to the number of fcf-disordered rungs in the 
ladder. The HERWIG event generator Q, as weh 
as the matrix element plus parton shower (MEPS) 
option of LEPTO are implementations of DGLAP 
evolution, with different approximate treatments of sub- 
leading terms. The ARIADNE generator |30t] uses the 
colour dipole model (CDM) in which disordered 

transverse momentum configurations are permitted but 
the radiated gluons are regarded as emitted by the 
proton remnant, so the connection with the picture in 
figure is not so clear. 

At very small classification of contributions 

according to the degree of ordering in x rather than 
kt becomes more appropriate, because logarithms of 
X dominate over those of Q^. The most widely used 
approach is that based on the BFKL equation p^ , 
which describes the a;-dependence of the structure 
function in terms of a sum over ladder graphs that 
exhibit strong ordering in x without ordering in 
transverse momentum. The original BFKL derivation 
was based on multi-Regge dynamics, in which kt values 
are limited. More recently, the same equation for 
the x-dependence of the structure function has been 
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Zt) of an emitted gluon, qt = qu/{l - Zi), ku = Yl]=i Itj 
is the corresponding exchanged transverse momentum, 
and ots = Cy^as/TT. The function l!^^{zi,kti) is a Regge 
form factor, 



/^n(zi,kti) = exp -as 



1 dz 
z 



dkt 



The CCFM expression corresponding to eq. 



(r) 



n 



d'^Qt dzi 



ZiasA^{zi,ku,qi) 



xQ{qi - Zi^iqi^i)0{qi - fi) 



(24) 
I is 
) 

(25) 



The main differences are the ordering condition qt > 
Zi^iqi-i, which corresponds to angular ordering of gluon 
emission due to coherence, and the replacement of the 
Regge form factor (e3) by the 'non-Sudakov' form factor 



Figure 7. Ladder graph for DIS at small x. 



derived by Ciafaloni, Catani, Fiorani and Marchesini 
(CCFM) from an approach based more closely on 
QCD, in particular on the coherence properties of the 
theory at small x. Although the CCFM approach leads 
to the same equation as BFKL when summed over all 
final states, it gives a different picture of less inclusive 
quantities Q. 

In order to compare the results of the two 
approaches, it is convenient to perform a Mellin 
transformation, introducing the moment variable lj 
conjugate to x: 



F{u;) 



dx x'^Fix) 



(21) 



Then the objective at small x is to sum to all orders in 
as the terms that are most singular at lu — 0. In both 
approaches the structure function is expressed as a sum 
over ladder contributions as in figure |^: 



oo 

F(^,Q/m)=^fW(u;,Q/m) 

=0 



(22) 



where r is the number of rungs in the ladder and /j, is 
an infrared cutoff. In the BFKL case one may write the 
ladder contributions in the form 



F, 



(r) 



n 

=f^ i=i 



d'^qt dz^ ^ 



z^asAj,{zi,ku) 



(23) 



where the variable qi is related to the transverse 
momentum qu and longitudinal momentum fraction (1 — 



An(zj, kti^qi) = exp -as 



dz 
z 



(26) 



This leads to the same asymptotic behaviour for the 
structure function: 



F„ 



Fn 



1 + 2as 



t 



2ai 



(27) 



where t = hi{Q/fi). However, the final state properties 
are different. For example, the mean number of ladder 
rungs, which gives the mean number of primary emitted 
gluon jets, is related to the quantity ^ rF^^"' , for which 
one finds Isl 



/ ^ ' BFK 



t 



2as 1- Gai — 



Y^rF. 



(r) 

CCFM 



2a 



s- +4a2 

UJ 



(28) 



Notice that the BFKL series for this quantity, like that 
for the structure function itself, involves at most one 
extra factor of 1/lu for each power of as, whereas the 
CCFM series contains up to two such factors. The 
new terms are subleading with respect to the number 
of logarithms of Q but dominate at small a: (w — > 0). 
The corresponding expansions for the mean number of 
primary emitted jets are 



{r)c 



as 




Q 










In 




1 - 


f as In 




X 
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Q 


14 


- as In^ 
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X 










X 



1-21 4 1 

-a„ in — 

6 " X 



(29) 



We see that more emission at small x is predicted by 
the CCFM treatment. The reason is that the angular 
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ordering prescription in eq. ( p5[ ) corresponds to an 
increased phase space compared to that available in 
eq. ( p3| ) for any fixed value of the infrared cutoff /i. 

So far, the CCFM approach has not reached the 
stage of detailed phcnomcnological application, and so 
the quantitative differences between the CCFM and 
BFKL predictions for final states in the accessible 
domain of x and are not yet clear. Predictions of 
CCFM and DGLAP evolution were compared in ref. 
at large down to x ~ 10~^ and the differences were 
not great; one would expect CCFM and BFKL to be 
even more similar. 

3.2. Transverse energy flow 

The new basic feature of small- a; physics (in both 
the BFKL and CCFM formulations) is the loss of 
strong ordering of exchanged and emitted transverse 
momenta. This suggests that one should see extra 
transverse energy emission in small-x events, especially 
in the region well separated from both the current 
and remnant jets. The HERA experiments do see a 
substantial transverse energy per unit rapidity in this 
region (about 2 GeV per unit rapidity) |^ |8), 
with some indication of a rise at small shown 
for example by the Hf data in figure |^. The curves 
show that standard DGLAP evolution as implemented 
in the MEPS model (defined above) does not produce 
such a rise, while the colour dipole model and a BFKL- 
based calculation |39) do, presumably because of their 
inclusion of /ct-disordered contributions. 

One point to bear in mind when comparing 
theoretical predictions for transverse energy fiow is that 
they are highly sensitive to the correct treatment of 
the hard scattering process that takes place at the top 
of the ladder in figure 0. This is because the mean 
emitted transverse energy is obtained by inserting a 
factor of kti into the integration in eq. ( p3|) or (|2^). 
As a result, the integral becomes more sensitive to the 
region ku ^ Q. In the case of the DGLAP evolution 
equation, which assumes strong ordering in kt, the kf- 
ordering above the i-th rung is in fact destroyed by 
inserting a factor of ku, and so the relevance of the 
equation becomes questionable. Preferably, one should 
compute the transverse energy emission at a particular 
rapidity from the full hard scattering matrix elements 
for producing additional partons in that direction. 

An illustration of the above point was provided by 
the HERWIG (version 5.7 and earlier) prediction for 
the transverse energy flow, which had a deficit in the 
rapidity range of interest. This was because the parton 
shower approximation to the hard matrix elements (with 
the choice of evolution variables used in HERWIG) left a 
region of phase space unpopulated by emitted partons. 
In version 5.8 a matrix element correction was included 
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Figure 8. Mean transverse energy flow per unit rapidity (at 
zero rapidity in the hadronic cm. frame), as a function of 
Bjorken x, for various ranges of . 



and the transverse energy flow became consistent with 
the data, without any retuning of parameters . 

3.3. Associated forward jets 

The most reliable indicator of new physics at small 
X still seems to be the associated production of 
jets in the forward (proton) direction with relatively 
large momentum fractions xjot and transverse momenta 
'Ztjot ^ Q [0 • In this region there is no phase space for 
DGLAP evolution with transverse momentum ordering, 
whereas the scope for BFKL evolution in x remains 
large. Thus a strong increase in production of such 
jets is expected with decreasing Bjorken x Since 
the cross section is fully inclusive with respect to what 
lies between the forward and current jets, the CCFM 
approach leads to the same conclusion. 

The HI data (figure H) |Q do indeed show a high 
event rate at small x. The data shown are for events 
at Q2 _ 20 GeV^ with jets with xjet > 0.025, 0.5 < 
g^jpt/Q^ < 4 and qt^ct > 5 GeV. The predictions of the 
colour dipole and MEPS models are also shown. While 
the colour dipole model does show a rise at small x, it 
predicts an even lower rate than MEPS if the cut on Xjct 
is increased to 0.05 \ ^T\ . 

At present the lack of a full BFKL-based event 
generator prevents a comparison with the BFKL 
prediction including selection cuts and detector and 
hadronization corrections. However, a parton-level 
BFKL calculation |E2| predicts a ratio of about 1.6 
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Figure 9. Forward jet production rate as a function of Bjorken 

X. 
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Figure 10. Comparison between conventional and 
instanton-induced DIS events. 



between the two x bins in figure ^, which is consistent 
with the data. 

4. Exotic QCD 

4-.1. Instanton-induced processes 

An intriguing possibiHty raised in recent theoretical 
papers |4^, but not so far subjected to 

experimental investigation, is that fundamentally non- 
perturbative features of QCD may lead to exotic 
new processes in DIS. Such processes can arise from 
the structure of the QCD vacuum, which contains 
topologically non-trivial gauge field configurations 
characterized by an integer value of a quantity called 
the Chern-Simons number, iVcs. Configurations with 
adjacent values of Ncs are separated by an energy 
barrier. The process of tunnelling through this barrier 
is represented by a space-time configuration called an 
instanton p6{ . 

An important feature of the instanton tunnelling 
process is that the change in the gauge field induces 
a change in the fermion sector of the theory. In the case 
of QCD the ferniionic change involves the axial charge 
Q5: 

AQ5 = 2AiVcs = 27V/ (30) 

where Nf is the number of active flavours. Since the 
axial charge measures the chirality or handedness of 
the ferniions, this implies that the process is chirality- 
violating. In electroweak interactions the effect is more 



spectacular: the tunnelling process leads to violation of 
baryon- and Icpton- number conservation ]4^ . Thus it 
is possible (though at present considered unlikely) that 
instanton-induced processes could be responsible for the 
net baryon number of the universe. This makes it all 
the more interesting to look for the analogous processes 
in QCD, which turn out to have much larger predicted 
cross sections, essential because as ~ 1/8 is much larger 
than tttu ~ 1/30. (For a tunnelling process, a factor of 
1/a appears in the exponent of the transition rate - see 
eq. dl) below.) 

Taking Nf —i, eq. (^^ implies that processes like 

7* ulUr , UR + g ULULdLSLULdhSL + gluons 

(31) 

can occur (figure |l^. The subscripts indicate that 
the outgoing quarks and antiquarks are all left- 
handed. The chirality violation would be difhcult 
to observe experimentally, but the relatively high 
multiplicity, including excess strange (possibly even 
charmed) particle production, should be characteristic. 
The typical number of gluons produced is expected to be 
of order 7r/2as ~ 10, further enhancing the multiplicity. 
Another distinctive feature is the relative isotropy of the 
final state in the rest-frame of the instanton subprocess, 
which transforms into a band about two units wide in 
the lab rapidity distribution, in contrast to the more 
normal jet-like configurations, as indicated in the lower 
half of figure |l^. 

A Monte Carlo generator for such events is under 
development B^: a typical simulated event is depicted 
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in figure |T^. The isolated particle at T^iab — —2.5 
is the electron, while the current-quark jet is around 
ryiab — —0.5. The densely populated band at ryiab — 2 — 4 
is the final state from the instanton-induced process. 

A general difficulty with instanton-induced processes 
is that they are normally suppressed by large inverse 
powers of Q'^ . The cross section is of the form 

.(.',Q'^)^lexp{--^.(.')} (32) 

where s', x', Q' are the cm. energy squared and Bjorken 
variables of the instanton subprocess (the blob I in 
figure |l^). S{x') is an increasing function, approaching 
unity at large x' but unfortunately not yet calculable 
at small x' . This translates into a contribution to 
the structure function F2 as shown in figure ^ The 
dashed curves are contours of constant S. The region 
labelled "data" corresponds roughly to the trend of the 
experimental data on F2. Thus at x < 0.3 a detectable 
signal could be present, if S becomes small enough. 



jet production rates; we can expect similar results from 
analyses of event shapes and jet fragmentation, plus a 
host of new data on power-suppressed corrections, which 
will improve understanding of the interface between 
perturbative and non-perturbative QCD. In the small- 
X region, theoretical ideas and experimental data on 
final state properties could be crucial in elucidating 
the mechanism that causes the rise in the structure 
function F2. A search for exotic final states due 
to instanton-induced processes could shed new light 
on non-perturbative QCD, and might establish the 
credibility of analogous electroweak processes as a source 
of baryon- and lepton-number violation within the 
Standard Model. 
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5. Conclusions 



Both experimental and theoretical studies of DIS final 
states are still in their infancy. Although things are 
not so simple as in e^e~ annihilation, because of 
the presence of the incoming proton and its outgoing 
remnant, DIS has the advantage of covering a wide 
range of at a single beam energy. This is already 
leading to evidence for the running of as from multi- 
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